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We have studied the effects of an external sinusoidal force in protein folding kinetics. The exter- 
nally applied force field acts on the each amino acid residues of polypeptide chains. Our simulation 
results show that mean protein folding time first increases with driving frequency and then decreases 
passing through a maximum. With further increase of the driving frequency the mean folding time 
starts increasing as the noise-induced hoping event (from the denatured state to the native state) 
begins to experience many oscillations over the mean barrier crossing time period. Thus unlike 
one-dimensional barrier crossing problems, the external oscillating force field induces both stabi- 
lization or destabilization of the denatured state of a protein. We have also studied the parametric 
dependence of the folding dynamics on temperature, viscosity, non-Markovian character of bath in 
presence of the external field. 

PACS numbers: 



I. INTRODUCTION 

Enhancement of reaction kinetics due to interplay be- 
tween barrier fluctuation rate and thermal noise-assisted 
barrier crossing events known as resonant activation^, 
is an interesting observation in early 1990s. This phe- 
nomenon provides a better understanding the mecha- 
nisms of various chemical and biological processes. Ex- 
amples include: dissociation kinetics of large molecules in 
coupled chemical systems^, oxygen binding mechanisms 
to hemoglobin 3 , modelling the dynamics of dye laser, 
ratchet models for the directional movements of molec- 
ular motors, transport through artificial nanopores^ etc. 
The interesting mechanism for enhancing rate via res- 
onant activation has triggered a numerous theoretical 
investigations for Markovian and non-Markovian nature 
of the bath and external deriving forces'*" 1 -^. This phe- 
nomenon has been experimentally realized^ in a tun- 
nel diode biased in a strongly asymmetric bistable state 
in the presence of two independent sources of electronic 
noise. The interference effects of resonant activation and 
stochastic resonance have also been studied to the aim for 
stochastic localization of particles confine d in a bistable 
potential as well as in a multi-well systemP^EH. 

The overwhelm majority of the previous studies on bar- 
rier crossing dynamics over a fluctuating energy barrier 
are based on the dynamics of a molecule having a single 
degree of freedom. It is modeled by a Brownain particle 
in a bistable or metatable potential with fluctuating bar- 
rier. But dynamics of the molecules having large number 
of degrees of freedom offers a significantly different situa- 
tion, due to its structural rigidity and several interaction 
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energies. Here, in addition to the energetic barrier and 
fluctuation statistics, the dynamics of the molecules is 
largely controlled by entropic factors. 

The folding dynamics of macromolecules like proteins 
is an example of thermally activated barrier crossing dy- 
namics in a multidimensional space. Here, all the degrees 
of freedom are coupled to each other through various in- 
teractions, such as, nearest-neighbor interactions, dihe- 
dral potentials, hydrogen bonds, ion pairs, van der Waals 
interactions etc. Over the years a considerable atten- 
tion has been focused on better understanding of folding 
mechanisms and potential energy landscapes. The stud- 
ies of protein folding dynamics under different external 
conditions such as salt concentration, temperature, con- 
finement, pH and viscosity of the medium et c 14 " 24 ! pro- 
vide important information about its structure and func- 
tionality. Again, the single molecule pulling experiments 
by highly sensitive force probes such as atomic force 
microscopy^H^ZI anc j optical and magnetic tweezers^ 28 * 30 l 
make it possible to realize various interactions due to in- 
ternal degrees of fre edom of a macromolecule. Motivated 
by the recent studies ^ * 26 * 28 * 30 ! in this context we have in- 
vestigated the effect of an applied external periodic field 
of non-thermal origin on the protein folding dynamics. 
We have assumed that the external electric field interacts 
with charge or dipole moment of the amino acid residues. 

To accomplish our goal we have co nsider ed the well- 
known coarse-grained off- lattice models^ 3 -^ for polypep- 
tide chain (where each amino acid residue is considered as 
a single bead centered at their C a position). The inter- 
action energies which play roles in the folding dynamics 
are taken into account by a Go-like HamiltoniarPsl. By 
setting the Langevin equations for each bead we have 
followed the dynamics of protein chains in presence of 
external fluctuations (field). 

A number of experimental studies^' 3 - 3 -' imply that 
Markovian dynamics cannot accurately account the effect 
of viscosity on the barrier crossing phenomenon in a so- 
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lution phase and the theory based on the non-Markovian 
dynamics shows a better agreement between theoretical 
and experimental results. Therefore, to capture the im- 
portant effects of non-Markovian dynamics and also for 
a sake of generality we have considered exponentially de- 
caying memory of the thermal fluctuations. Moreover, 
the present study is an extension of earlier studieipSEID 
on barrier crossing dynamics in low dimensional systems 
in presence of periodic force. To compare the features 
of folding kinetics of macromolecules with one dimen- 
sional barrier crossing problems, we also simulate a non- 
Markovian Langevin dynamics in a bistable potential in 
presence of barrier fluctuations. 

Specifically our objective here is twofold. First we in- 
tend to explore the effects of an applied external oscil- 
lating field on protein folding to the goal of extracting 
generic features of multidimensional barrier crossing dy- 
namics in contrast to the one dimensional cases. The sec- 
ond objective is to investigate the effects of the oscillating 
time periodic force field on the following three important 
features of protein folding kinetics: (i) turnover behavior 
of mean folding time with solvent viscosity, (ii) U-shaped 
mean folding time vs. temperature profiles, and (iii) dou- 
ble minimum in mean folding time as a function of the 
correction of thermal fluctuations. 

To address the above challenging issues we simulate 
folding dynamics of the 16-rcsidue peptide /3-hairpin (C- 
terminal from protein G, PDB ID: 2gbl) and 76-residue 
protein ubiquitin (PDB ID: lubq). The PDB structures 
of these proteins are shown in Fig. 1. 
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FIG. 1: (Color online) The PDB structures of two proteins 
studied in this work. For the cutoff distance d c — 6.5 A , 
the total number of native contact is equal Q ma x = 13 and 
99 for hairpin and ubiquitin , respectively. 



tractive and non-native interactions are repulsive. The 
energy of a configuration of a protein is specified by the 
coordinates of the C a atoms and is given by^ 



E = ^ K r (r iti+ 
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r i.i+i is the- distance between i-th and (i + l)-th beads. 
9i is the bond angle formed by three subsequent beads: 
(i — l)-th, i-th and (i + l)-th. 0i denotes the dihedral 
angle around the i-th bond and is the distance be- 
tween the i-th and j-th residues. The subscripts 0, NC, 
and NNC refer to the native configuration, the native 
contact, and the non- native contact, respectively, r^ij is 
the distance between the i-th and the j-th residues in the 



native conformation and = 



Amino acid residues 



are assumed to be in the native contact if r^j is less than 
a given cutoff distance(d c ). Here, we have assumed the 
cutoff distance d c = 6A. 

The first term of Eq.Q presents harmonic potential 
due to chain connectivity between two adjust beads. The 
second term is also harmonic potential arising due to 
bond angle between three subsequent beads. The third 
term presents dihedral potential for every four adjacent 
C a atoms. Dihedral potential is a sum of two peri- 
odic components with periods: = 2ir/((f) i — (f) 0i ) and 
T0 = 27r/3(0i — c/)Qi). The last two terms in Eq.(fT|) 
are due to the non-local native interactions and the 
short-range repulsive force for non-native pairs, respec- 
tively. Parameters K r , Kg, K^, en denote the relative 
strength of each kind of interaction; we choose K r — 



lOOen/1', Kg = 20e H /rad 2 , K^ = e H , K^> = 0.5e ff , 
where en is the characteristic hydrogen bond energy and 
C = 4l. 

The dynamics of the protein chain in a thermal bath 
can be described by setting a Langevin equation for each 
bead ( C Q -atoms of each amino acid residue) of the pro- 
tein chain. 



,(3) 



II. THE MODEL 

To study protein folding dynamics we consider coarse- 
grained off-lattice models^ for polypeptide chains in 
which each amino acid residue is represented as a 
single bead centered at its C a position. Moreover, 
we follow the d yna mics of polypeptide chains by Go- 
like HamiltoniarPS, in which the interactions between 
residues forming native contacts are assumed to be at- 



-y(t - t')rd£ + rf{t) 



(2) 



where m denotes the mass of a bead and F c — VE cor- 
responds to force derived from the potential energy asso- 
ciated with the protein molecule (Eq. ([!])). The potential 
force F c depends on the positions of the all beads. There- 
fore, TV coupled Langevin equations are needed to fol- 
low the time evolution of a protein having TV amino acid 
residues (beads). The thermal fluctuations due to en- 
vironment are modelled by colored Gaussian noise. The 
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frictional kernel 7 (t) is related to thermal fluctuation £(t) 
by the well-known fluctuation-dissipation relation. 

(Vk(t)m(t')) = k B T6 kn (t-t') (3) 

where, k B and T are the Boltzmann constant and tem- 
perature of the system, respectively, rjk is the fc-th com- 
ponents of ff(t), where k, I = x, y, z\ the three orthog- 
onal components. To capture essential features of the 
non-Markovian dynamics we consider exponentially de- 
caying frictional memory kernepsHID w ith the following 
form: 

7 (t-t') = ^exp(-^) (4) 

T T 

where 70 is the frictional coefficient in the Markovian 
limit and r bears the memory effect of the non-Markovian 
dynamics. Then r)(t) is a solution of the following differ- 
ential equation: 

i] = -r]/T-\ CO)- 

T 

Where, £(t) is a Gaussian white noise having variance 
two. It should be noted that for the frictional memory 
kernel (JI]) the integro-differential ([2| can be simplified as 

mf= F c + rf(t) and rf = -ff/r - j f/r + ^ l0 * B ' r ({t). 

We have considered the situation where the dynamics 
polypeptide chains is affected by an external sinusoidal 
force. It implies that free energy of the system periodi- 
cally oscillates in an asymmetric way. Under such situa- 
tion the dynamics of the protein chain is governed by the 
following equations 

mf=F c - [ 7 (t - t'ydt' + ff(t) + A sin ut (5) 
Jo 

Unlike the standard pulling experiments (where the ex- 
ternal force is applied to termini of bio- molecules), here 
the external oscillatory force filed acts on all amino acid 
residues. However, the last term in the above equa- 
tion accounts effective interaction between the charge or 
dipole on the amino acid and electric field. The resultant 
force vector corresponding to this interaction affecting 
the acceleration vector is parallel to the r. Thus A in 
the above equation is parallel to the acceleration vector. 
We have implemented this in the simulation scheme. 

A. Simulation 

In order to analyze the effects of an external field on the 
protein folding dynamics we have calculated the folding 
time for different parameter sets (frequency of the exter- 
nal field, dissipation constant of the medium, tempera- 
ture, noise correlation time etc.) by numerically solving 
the Langevin equation^ . The relevant equations^ and 
([5]) were integrated using the velocity form of Verlet al- 
gorithm with the time step At = O.OOlri, T£ is the char- 
acteristic time scale of the system which are defined by 



t~l = (ma 1 1 en) 1 ! 2 w 3 ps, where a is the characteristic 
bond length between two successive beads a = 4A. The 
mean folding time is the averaged first passage time from 
a random configuration (thermodynamically unstable) to 
the native state (thermodynamically stable state). Our 
calculated mean folding time is the averaged over 250- 
1000 trajectories depending on the values of different pa- 
rameters to have good statistics. To this end we would 
like mention that in the present paper we have exploited 
the numerical scheme of our earlier study^. 

B. Results and discussions 

We have calculated the mean folding time(ty) of 
/3— hairpin as a function of frequency of the external 
field at temperature T — 0.53e#/fcs = 298K, where 
e.H = 0.98 K cal/mol is the hydrogen bond energy. This 
is depicted in Fig. 2. The folding time of 0— hairpin first 
increases with increasing frequency of the external force 
followed by its decrease after passing through a maxi- 
mum. With further increase of frequency the folding time 
starts increasing. Thus, the mean folding time versus 
driving frequency profiles possess both a maximum and 
a minimum. A similar feature has been observed (not 
shown here) at the Markovian limit, r — > 0. It should 
be noted that the maximum is not observed in the one- 
dimensional barrier crossing problem. To demonstrate 
this issue and for a qualitative comparison with an one- 
dimensional case, we have considered barrier crossing dy- 
namics of a Brownian particle, modelled by the following 
generalized generalized Langevin equation 4 ^, 

mv = q-q 3 - [ 7 (t - t')v(t')dt' + £(t) (6) 

where q and v are the position and velocity of a Brownian 
particle. Here, the folding time is defined as a mean bar- 
rier crossing time, is obtained by solving the above equa- 
tion and plotted in Fig. 3. The results for the correspond- 
ing Markovian limit have been depicted in the inset of 
Fig. 3. Fig. 3 reveals that, for an one-dimensional case, the 
mean escape time versus driving frequency plots possess 
no maximum. This aspect has also been reported ear- 
lier by several groups using numerical experimenlP^ and 
analytical calculation^. Also, there are experim ental 
and theoretical studies modelling catalytic reaction d 35 ! 44 ! 
which show that mean escape time vs. driving frequency 
plots exhibit only a minimum but no maximum. Thus, 
appearance of the maximum at low frequency regime is 
a generic signature of barrier crossing dynamics in many 
dimensional systems. In protein folding kinetics, in ad- 
dition to energetic barrier, entropy of activation plays 
a significant role. The appearance of the maximum (in 
the tf vs. lo) at low frequency regime may be attributed 
to increase of the fluctuations in configurational entropy 
during the barrier crossing rather than activated trans- 
port towards the transition state by the periodic force. 
For further increase of frequency the energy transfer to 
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FIG. 2: (Color online) Plot of mean folding time tf of /3— 
hairpin vs the frequency of the external driving force. We 
chose T = 298A", A = 0.3 pN, r = 1.0 ps. The inset presents 
the results of Markovian limit, r — > 0. 



activate towards the transition state dominates over the 
increase of fluctuations in configurational entropy and 
the mean folding time decreases until the energy transfer 
rate becomes small. Thus the maximum and the min- 
imum in tf vs w plots is a result of interplay between 
two important quantities, entropy and activation energy. 
Similar features of folding time as a function of driving 
frequency has been observed for the higher driving am- 
plitude and damping. It has been presented in Fig. 4. 

Figures 2 and 4 reveal that the position of the min- 
imum in the mean folding time vs. frequency plots is 
very sensitive to viscosity (damping) of the medium. Po- 
sitions of the minima are shifted to the lower frequency 
with increasing damping. The pleasurable explanation 
of this behavior is as follows: At low viscosity, because 
of the strong spatial diffusion the acceleration of folding 
kinetics starts to work over the entropy effect slowly com- 
pared to high damping case. Therefore, the minimum at 
the higher driving frequency for lower damping constant 
is observed. 

Our another interesting observation from the Figs. 2 
and Fig. 4, the minimum is shallow and flat for low damp- 
ing cases compared to the high damping situation. Be- 
cause of strong spatial diffusion in the former case, the 
effect of the periodic force (by virtue of doing mechanical 
work on the particle in the dynamics) is weaken and it 
suppresses activation as well as slow down the variation of 
folding time with increase of the driving frequency. Thus, 
the increase of the damping constant enhances activation 
at the cost of decreasing its robustness. 

To generalize our conclusions we have also studied fold- 
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FIG. 3: (Color online) Plot of mean first passage time MFPT 
vs ui the frequency of the external driving force. We chose 
the double well potential as V(q) = l/4g 4 — l/2q 2 . We use 
T = 0.1, Ao = 0.25, r = 1.0.. In the inset same plot is drawn 
in the limit r — > 0. (units are arbitrary) 
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FIG. 4: (Color online) This plot shows the dependence of 
mean folding time of ft— hairpin on the frequency of the ex- 
ternal derive in the heavy damping situation. The chosen 
parameters are T = 298A", r = 1.0 ps, Ao = 1 pN. 

ing kinetics of a long protein, ubiquitin. Figure 5 presents 
mean folding time of ubiquitin versus driving frequency 
for different values of the damping constant. Here also 
the mean folding time passes through a maximum and a 
minimum. Other behaviors are also very similar to /3— 
hairpin. 

Next, we have explored effects of the non-Mar kovian 
bath in the folding dynamics in presence of external field. 
The folding time of /3— hairpin shows two minima with 
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FIG. 5: (Color online) Plot of mean folding time tf of ubiq- 
uitin vs the frequency of the external driving force. We chose 
T = 298K, A = 0.1 pN, r = 1.0 ps. 
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FIG. 7: (Color online) This plot shows the dependence of 
mean folding time of /3— hairpin on the viscosity of the 
medium. The chosen parameters are T = 298-RT, r = 
1.0 ps, A = 1 pN. 
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FIG. 6: (Color online) This plot shows the dependence of 

mean folding time of (3— hairpin on the correlation time of FIG. 8: (Color online) This plot shows the dependence of 
the non-Markovian noise. The chosen parameters are T = mean folding time of j3— hairpin on the Temperature of 
298i£", 7 = l.Ops -1 , Aq — 1.0 pN. the system. The chosen parameters are 7 = 2 ps -1 , r = 

1.0 ps, A Q = 1 pN. 



noise correlation time (shown in Fig. 6). The reason for 
appearance of unusual two minima has have discussed in 
detail in our recent studj3^. However, the second min- 
imum appears at higher correlation time in presence of 
the periodic force. To explain this we recall the fluctu- 
ation dissipation relation ([3]). It implies that the vari- 
ance of noise decreases with increase of noise correlation 
time. Increase of noise correlation time leads to decrease 
of fluctuations in entropy. Because of higher entropy in 
presence of the driving force, the noise correlation starts 
playing role in the dynamics at its larger value. Thus, 
the shifting of positions of the minima are attributed to 
the excess entropy due to the external field. 



Figure 7 depicts the variation of folding time of 
/3— hairpin as a function of viscosity of the medium in 
presence of external periodic force. We observe that the 
folding time differs significantly around only the region 
where the folding kinetics is accelerated due to the ap- 
plied periodic field. 

To elucidate the effect of external forcing on the U- 
shaped mean folding time versus temperature profile, we 
calculate the folding time both in presence and absence 
of the external force. This is depicted in Fig. 8. It shows 
that the folding time is greater than the corresponding 
unperturbed case in the regime where increase of entropy 
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dominates over the energy transfer to activate towards 
the transition state and is lower when the latter dom- 
inates over the former. Acceleration of folding kinetics 
due to the external field causes a shift of the minimum 
(see Fig. 8) towards the higher temperature and also low- 
ers the width of U-shaped temperature profiles. 

III. CONCLUSION 

Based on the coarse-grained off-lattice models for a 
polypeptide chain and setting Go-like Hamiltonian for 
the system we have followed the dynamics of the protein 
chain by solving N— coupled Langevin equations in pres- 
ence of an external field. In order to make the model 
more realistic we have assumed the non-Markovian heat 
bath. Our main conclusions of this study are as follows: 
(A) Under influence of an oscillating electric filed the 
mean folding time of /?— hairpin shows a maximum and 
a minimum in mean folding time vs driving frequency 
plots. This feature is in a sharp contrast to the one- 
dimensional barrier crossing problem where only a min- 
imum is observed in the same profile. Thus, periodic 
force can induce stabilization or destabilization of the 
denatured state of a protein. The folding kinetics of a 
longer protein, ubiquitin, also exhibits a similar feature. 
Therefore, the above observation is true for both long 
and small proteins. 

(B) Even in presence of an oscillating force field of am- 
plitude 0.1-1 pN the following three important features 



of protein folding kinetics remain intact: (i) turnover be- 
havior of mean folding time with solvent viscosity, (ii) U- 
shaped mean folding time vs. temperature profiles, and 
(iii) double minimum in mean folding time as a function 
of the correction of thermal fluctuations. It implies that 
the above properties is very robust to external perturba- 
tions. But the following new features are noticed for the 
presence of the external force field. 

(a) In the presence of the periodic force the second min- 
imum appears at relatively higher noise correlation time 
of the thermal noise in the plot tf it vs. r. This is due 
to excess entropy of the system, which is produced by 
external fluctuations. 

(b) Increase in the damping strength enhances accelera- 
tion of folding kinetics by the periodic force at the cost 
of decreasing its robustness. 

(c) Finally, we observe acceleration of the folding kinetic 
due to presence of the external field in the variation of 
mean folding time with viscosity of the medium and tem- 
perature. 

We hope our theoretical findings could be verified ex- 
perimentally measuring protein folding time in the pres- 
ence of a monochromatic isotropic electromagnetic radi- 
ation. For a given intensity of the radiation amino acid 
residues experience a resultant amplitude vector Aq. In 
presence of an isotropic electric field the magnitude of A a 
would be independent on the orientation of the protein 
chain and time. We have considered this aspect in our 
present study. 
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